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Abstract 

We discuss the threshold anomaly in ultra-high energy cosmic ray physics 
by assuming that the matter world just be perturbation of a de Sitter universe, 
which is consistent with the recent astronomical observations: about two thirds of 
the whole energy in the universe is contributed by a small positive cosmological 
constant. One-particle states are presented explicitly. It is noticed that the 
dispersion relation of free particles is dependent on the degrees of freedom of 
angular momentum. This fact can be regarded as the effects of the cosmological 
constant on kinematics of particles. 

PACS numbers: 98.70.Sa, 95.30.Cq, 95.85.-e. 



1 Introduction 

The origin of the ultra high energy cosmic rays (UHECRs) is one of the outstanding 
puzzles of modern astrophysics. Present understanding of the phenomena responsible 
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for the production of UHECR is still limited. Currently, there are generally two sce- 
narios to produce the UHECRs. One is the "bottom-up" acceleration scenario with 
some astrophysical objects as sources^|2]. The other is called "top-down" scenario 
in which UHECR particles are from the decay of certain sufficiently massive particles 
originating in the early Universe [KI • 

Decades ago, Greisen, Zatsepin and Kuz'min (GZK) j3] discussed the propagation 
of UHECRs through the cosmic microwave background radiation (CMBR). Due to pho- 
topion production process by the CMBR, the UHECR particles will lose their energies 
drastically down to a theoretical threshold, which is about 5 x lO^^eV. The mean free 
path for this process is only a few Mpc [Hj. This is the so-called GZK cutoff. However, 
we have observed indeed hundreds of events with energies above lO^^eV and about 
20 events above lO^'^eV jOl-fn]- The puzzle can be considered to be some cosmic ray 
threshold anomaly: energy of an expected threshold is reached but the threshold has 
not been observed yet. Although some physicists argued that the GZK cutoff might 
actually have already been seen it need more evidences and careful analysis to 
become the consensus in scientific circles. In this paper, we will take the existence of 
UHECR paradox as an acceptable working assumption, on the basis of which we try 
to unravel the puzzle. 

The cosmic-ray paradox was suggested to be solved by the departures from Lorentz 
symmetry in the end of the last century fT^ IT^ . The basic idea is to introduce tiny 
Lorentz-violating terms into the standard model Lagrangian and derive a parameter- 
dependent threshold. As long as the parameters are taken some appropriate values, the 
threshold will be above the energies of all observed events. A related but alternative 
scheme is proposed ^1], in which Lorentz symmetry is broken by Planck-scale effects. 
Another approach is to assume deformed Lorentz symmetry, constructing the theory of 
doubly special relativity (DSR) [13 CEl IHj • In the approach, a new dispersion relation 
is written down and then an enhanced threshold is obtained. 

The recent astronomical observations on supernovae ^HJ and CMBR [20] show 
that about two thirds of the whole energy in the Universe is contributed by dark energy. 
The simplest model for dark energy is that it is devoted by a small positive cosmological 
constant (A). Then, the Universe can be regarded as a de Sitter (dS) spacetime in the 
zeroth approximation. The physics in dS spacetime has been discussed extensively 
(see, for example, |2T|-[23]). 

In our previous paper [21j, we have discussed the possibility of the cosmological 
constant as the origin of threshold anomaly in a dS spacetime. In a kind of simplified 
case, we obtained a positive conclusion. However there are two imperfects in Ref. 
which are also motivations for this paper. Firstly, although the dominated ingredient in 
the Universe is dark energy, the matter also contributes about 1/3 to the density of the 
Universe. The challenging task is to find a way of dealing with the universe dominated 
by the dark energy and matter. At least at the present stage of the evolution of the 
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universe, one has few reasons to omit effects of any of the two ingredients. Therefore, 
we need to investigate the Universe including matter in order to more approach to the 
Universe nowadays. Secondly, in |21] , only some simplified cases are considered, which 
are not necessary the true ones in our Universe and will be improved on in current 
paper. 

In this paper, we present a scenario in which the dark energy is supposed to come 
from the cosmological constant and that the spacetime is a dS one, and the matter in the 
Universe is dealt with as a perturbation around the dS background. The perturbation 
of the spacetime geometry is derived firstly. Then we discuss the motion of a free 
particle in the asymptotically de Sitter Universe, and its kinetics is set up primarily. 
Meanwhile, we derive a general form of dispersion relation for free particles moving 
in the Universe. This formalism is used to investigate the UHECR propagating in 
the cosmic microwave background. We obtain explicitly the corrections of the GZK 
threshold for the UHECR particles interacting with soft photons, which are dependent 
on the cosmological constant as supposed in the beginning of the paper. We show 
how the threshold varies with a positive cosmological constant and additional degrees 
of freedom of the angular momentums of interacting particles. It should be noticed 
that, for a positive cosmological constant, the theoretic threshold tends to be above the 
energies of all the observed events. Thus, we may conclude that the tiny but nonzero 
cosmological constant is a possible origin of the threshold anomaly of the UHECR. 

The paper is organized as follows. In Section 2, we discuss the Friedmann equation 
for a constant curvature spacetime with a homogeneous density perturbation. An 
explicit solution including metric and corresponding connection is presented. The 
section 3 is devoted to the investigation of kinematics. Approximate conservation laws 
of momentum and angular momentum are obtained along the geodesies. By solving 
equations of motion of a free particle, we obtain a remarkable dispersion relation, which 
includes degrees of freedom of angular momentum. In Section 4, by taking effects of a 
tiny but nonzero positive cosmological constant into account, we show that the theoretic 
threshold is above the energies of all the observed UHECR events. In the last section, 
we present conclusions and remarks. 

2 Perturbation of the Friedmann equations 

The Einstein equation with cosmological constant is 



- 7;9abR - ^9ab 



SttGT, 



ab 




The stress-energy tensor in Eq. can be of the general perfect fluid form 



Tab = (P + P)UaUb - PQab , 



(2) 
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where p, p are the proper pressure and energy density, respectively. 

While observations of the distribution of galaxies in our Universe show clustering of 
galaxies on a wide range of distance scales, on the largest scales the galaxy distribution 
appears to be homogeneous and isotropic, namely, p and p are only functions of variable 
t. The homogeneous and isotropic universe is described by the Robertson- Walker 
metric 



ds"^ = dt^ - a?{t) 



(3) 



1 — kr"^ 

where k denotes spatial curvature of the spacetime. 

Thus, the general evolution equations for homogenous, isotropic universe read 

3a = —AnGip + ?>p)d 

.. .2 (4) 

aa + 2a +2k = AtxG{p — p)aP , 

where p = Pa + Pm, Pa and pm denote the energy density contributed from cosmological 
constant and matter, respectively. 

The stress-energy tensor should be covariantly conserved, i.e., T"*^ = 0, which 
reduces to 

P = -3?(p + p). (5) 
a 

It should be noticed that pm/ Pa is a small parameter for the evolution of the Uni- 
verse (at least in the present stage). To discuss deviation of physical quantities in 
an asymptotically de Sitter spacetime from those in a dS one, we would introduce a 
small parameter e characterizing the effect of the matter on the spacetime geometry. 
Then, one can assume the scale factor a{t) has a perturbation by the matter density 
Pm around a dS background, 

a{t) = a{t) [1 + ef{t)] = R cosh(t/i?) [1 + e/(t)] , (6) 



where R := ^^/3/A. 

From the first approximation of the covariant conservation equation 

p + 3H{p + p)=0, (7) 

and 



one can obtain 



fit) = ^ (tanh(t/i?) cot-^[sinh(t/i?)] - ^^^^^) • (9) 
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In the derivation, the asymptotic condition hm f{t) — ^ has been used. The asymptot- 

t— >oo 

ically dS spacetime with the above form of matter can be reahzed as a four dimensional 
hypersurface embedded in a five dimensional fiat space 

ds' = {dir - {d^r - {dir - id^r - {dir , (lo) 

such that 

- - - - = -R\t) . (n) 

^'^ (/i = 0, 1, 2, 3, 4 ) is related to the coordinate (t, r, 9, (p ) in Eq. (jH]) as follows 

m 

^0 - eO(t = 0)= f ^l + d\T)dT , = + ef) , {A = l, 2, 3, 4) , (12) 

^0 



where 



C,^ = r cosh(t / R) sin 6 cos (f), 
= r cosh(t / R) sin 6 sin 0, 
= r cosh(t/i?) cos 9, 



(13) 



^4 _ ^R2 ^r^cosh{t/R). 
Applying Eqs. © and to the expression of in (fT^ . one obtains 

= e° + eR'^^^fit) + epoi?^ (coth(2t/i?) - 1) , (14) 

where 

e° = i?sinh(t/i?). (15) 

Similarly, in the course of deriving Eq. (fT^. the condition lim = is used. i?(t) 

t— ►oo 

in Eq. (fTT|) reads 

R{t) = R{1- epoR^ smh{t/R){coth{2t/R) - 1)) . (16) 

It is obvious that 

lim R{t) = R . (17) 

t— ►oo 

For convenience, we write the metric in terms of the Beltrami coordinates x"" = 
R-^ , (^^ 7^ , a = 0, 1, 2, 3) l^nj-EHl- Due to the additional matter, the metric 
of the 4-dimensional spacetime in Beltrami coordinates has a perturbed form 

ds^ = Qabdx^dx^ = {gab + ehab)dx''dx^ , (18) 
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(19) 



where gab is the metric in the empty dS spacetime with the form 

gab = \ — , Vab = aiag(l, - 1, - 1, - 1) , 

a := a{x, x) = 1 —5 — > 

and hab is the perturbation from the additional matter. 
From Eq. (fTSj) . the following results can be obtained 

hoa = , 

= -Vi^f-^^^ '9f, = -2f{t)R~^a\t) , j = 1, 2, 3 ; 
where ^gf^ is the Beltrami metric on 3-sphere 

3 ^ ^ij ^ik^jlX X 

= ^ " ' (21) 

0-3(0;) = 1 + R~'^6ijx''x^ . 

In order to investigate the motion of free particles, Christoffel coefficients are listed 
as follows 



(20) 



= ^tanh(t/i?)cot'i[sinh(t/i?)]— . 
2 0-3 



R^aix) 



T% = ^{cosh(t/i?)[l + tanh2(t/i?)]cot-i[sinh(t/i?)] -2tanh(t/i?)}(Ti/2 3^B 



- § tanh' {t/R) cot-^ [sinh(t/i?)] '''''' 



= ;^cot-^[sinh(t/i?)]5i, 
ff. = -^tanh(t//2)cot"^[sinh(t//?)] 



X 



^5^ + x^d'' 



2 ^ ' ' ' " as 

(22) 

It is not difficult to calculate the deviation of the geodesies in an asymptotically dS 
spacetime from those in the empty dS one 

= + = a- (-—] + eP'^ —— (23) 

ds^ '"^ ds ds ds\a ds I '"^ ds ds 
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In the asymptotically dS spacetime with homogeneous matter, one can define the 4- 
momentum formally as 



1 dx" 
a ds 



and E:=P'^ . 



From Eq. (j23|) . one has 



dP" 



-eVl^P^dx^ 



where v'^ is a constant defined as 



dx"' 



(24) 



dx° ■ 



For the particle moving along a one-dimensional curve, one can set 



X 





without loss of generality. That is, one can discuss the motion of a free particle in 
x^ — x^ plane. In this case, the deviations of energy and momentum in the perturbed 
dS universe from those in the empty dS universe can be calculated numerically. In 
Fig. 1, we present a plot for the relative errors of energy and momentum vs. the time 
coordinate x° (we use the convention c = 1 throughout the paper). 




100 



X (Mpc) 



Figure 1: The solid line describes the relative error of one-dimensional momentum for 
a free particle moving in the perturbed dS universe, while the dashed line describes its 
relative error of energy. 

These results show clearly that, as a feasible approximation, we can substitute 
energy and momentum in the empty dS universe for the counterparts in perturbed 
one provided that the distance between source of UHECRs and the earth is less than 
lOOMpc. 
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3 Kinematics around dS spacetime 

In the asymptotically dS spacetime with homogeneous matter one can still define 
the five dimensional angular momentum M^'^ of a free particle with mass mo as the 
form 

M- = ™„(ef -ff) . (25) 

where s is a parameter along the geodesic. In the Universe, there is no transla- 
tion invariance and so that one can not introduce a momentum vector. However, 
it should be noticed that, at least somehow, one may define formally a counterpart of 
the 4— momentum P for a free particle in the Universe as 

:= R-^M'"' = . (26) 

as 

In the same manner, the counterparts of the four dimensional angular momentum J"* 
can be assigned as J"'' := M°*. 

It is easy to show that ^^{:= a{x, x)~^^'^x^) is invariant under the spatial transfor- 
mations. Thus, we can say two spacelike events are simultaneous if they satisfy 

a{x, x)-ix° = ^° = const. . (27) 

Therefore, it is convenient to discuss physics of the asymptotically dS spacetime in the 
coordinate {C,^, x"). In this coordinate, the metric can be rewritten into the form 



Ap2)2 1 + Ap2 



28) 



3 

:2 



where p"^ '■= 2;"a;" and dQ^ denotes the metric on 2-dimensional sphere S"^ 

a=l 

The Klein-Gordon equation describes motions of a scalar field (f){x), 

(□ - ml) 0(x) = , (29) 
where □ is the d'Alembertian operator defined as 

□ := --L9,(v^rt) • (30) 



-9 

The Dirac equation for the spinor field \l/(x) has the form[^ 

H7"(5a-rj+mo]( ) = 0, (31) 

where is the Ricci rotational coefficient [22] • For free spin 1/2 particles, it is easy to 
check that the components tpa{x) of a spinor satisfy the relation 

{n - ml) ipaix) = , a = 1,2. (32) 
8 



In Lorentz gauge, one can simplify the Maxwell equation without source as 

UA" = . (33) 

From Eqs. (j^^ . (jSH) and ()33|) . one knows that scalar fields and components of spinor 
and vector fields can be described uniformly as (mo = for vector field) 



'-9 



(34) 



In the coordinates (.^°, x'), one can rewrite the d'Alembertian operator as the 
following form 

^5a(v^rt) = - (1 + A^r) 4 - (4Ae° + 3e(l + A^T)^) 

+(1 - 26/) (1 + xeey (i + vY k + sp^^^,] (35) 
H^-2ef)[i+xeer\i+v)p~'di, 

where A = -R~^, uu' = 1 and denotes the Laplacian operator on 5^. 

To solve the equation of motion, one writes the field $(.^'',x*) into the form 

$(e°,p,u) = T(e°)f/(p)i,^(u) . 

Thus, one transforms the equation of motion into OU] , 

;i + \eef{i + 26/)9|o + {^xeii + \ee){i + 26/) + 36^(i + xee?^ o^o 

+ mlil + 26/)(l + Xee) + ie' - ml)] T(e°) = 0, 



2 ^2 



:i + Ap2)2 p2(l + ^^2)J_ 



U{p) = 0, 



[dl + l{l + l)]YiM=0, 



(36) 



where Yim{u) is the spherical harmonic function and e is the constant from separating 
variables. 

For the irrelevance of the expression of r(^°) to our discussion, we can focus our 
attention on the last two equations. The solutions for the radial equation of the field 
is 

U{p) ~ p'(l + Xp')'/'F Q(/ + k + 1), ^-{l + k), 1 + ^; - Ap^) , (37) 
where k denotes the radial quantum number 



e - 2k - X'\e^ - ml) =0 
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To be normalizable, the hypergeometric function in the radial part of the wavefunction 
has to break off, leading to the quantum condition 

l + k = -2n, (2neN). (38) 
Then, one obtains the dispersion relation for a free particle moving in the Universe 

=ml + e''^ + \{2n + l){2n + l + 2) , (39) 
where the term e'"^ denotes what is independent of the parameters n and /. 



4 UHECR threshold 

In this section, we investigate the UHECR threshold in the covariant framework of 
kinematics in an asymptotically dS spacetime set up in the preceding sections. 

One considers the head-on collision between a soft photon of energy E^, momentum 
q and a high energy particle mi of energy Ei, momentum pi, which leads to the 
production of two particles m2, with energies £'2, E^ and momentums p2, Ps, 
respectively. From the energy and momentum conservation laws 

El + E^ = E2 + E3 , 

(40) 

Pi-q = P2+P3 ■ 

In the C. M. frame, m2 and are at rest when the threshold is reached, so they have 
the same velocity in the lab frame and there exists the following relation 

^ = ^ . (41) 

It is convenient to use the approximate formulae of dispersion relations for the 
soft photons and the ultra high energy particles 

E', = q' + \;, (42) 



E, = ^m^ + pl + \*c^p, + ^ + ^^, (^ = 1, 2, 3). (43) 

where A* := A(/^ + 2n^)(/^ + 2n^ + 2) ^ A/^(/^ + 2) and X* := X{li + 2ni){li + 2ni + 2) k. 
+ 2) with the conjecture that I ^ n . 
The obtained threshold can be expressed as the form 

(m2 + -ml + \*(l + ^)+ X; + XI 
' 2{E, + ./E^) ■ ^^^^ 



10 



The usual GZK threshold could be recovered when the parameter A*, which is depen- 
dent on the cosmological constant, runs to zero. 

The conservation law of the angular momentum imposes a constraint on the pa- 
rameters A*, 

At + a; + 2ALi ■ = a; + A* + 2AL2 ■ L3 . 

Making use of the relation, one can rewrite the A* dependent terms of the threshold as 
the following 

+ Kt^, +K + 2ALi ■ L, - 2XU ■ L3 

If A2 and A3 take value of the same order with A* (less than the square of energy of a 
soft photon), the A* dependent terms can be omitted j2lj. We will investigate the case 
of A2 + A3 > A* , and the threshold (jHj) is of the form 



{m, + msf - m? + A^^ + A*^ - 2XU ■ L 



th, A 



3 



(46) 



Now, we can study the photopion production processes of the UHECR interaction 
with the CMBR 

p + 7 ^ p + vr . 
The corresponding threshold for this process is given by 

^UHECR ^ (rn^ + m.Y -ml + A^^ + A;^ - 2AL^ • 
'''' - 2{E,+ ^E^) ■ ^^'^ 

To show the behavior of the threshold in the A*-parameter space clearly, we should 
discuss some limit cases in detail. 

In the case that the out-going nucleon has zero angular momentum, the threshold 
dHj) reduces as 

we provide a plot for the dependence of the threshold E^^^^ on the cosmological 
constant and angular momentums (the in-coming photon and out-going pion) in Fig. 
2. 

In the case that the out-going pion has zero angular momentum, the UHECR 
threshold takes the form 

We present a plot for the dependence of the threshold E^^^^ on the cosmological 
constant and angular momentums (the in-coming soft photon and out-going nucleon) 
in Fig. 3. 
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A;(xlO""Gev2) 10 



Figure 2: The cosmological constant and angular momentums (of in-coming soft photon 
and out-going pion) dependence of the threshold E'^^^f'^^ in the interaction between 
the UHECR protons and the CMBR photons (A* in units of m^/10) . 



A*(m2/10) 




A;(xlO""Gev2) 10 



Figure 3: The cosmological constant and angular momentums (of the in-coming soft 
photon and out-going nucleon) dependence of the threshold E^^^^ in the interaction 
between the UHECR protons and the CMBR photons (A^ in units of m^/10) . 
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A;(m2/10) 




Figure 4: The cosmological constant and angular momentums (of the out-going nucleon 
and pion) dependence of the threshold -E^^^f *^ in the interaction between the UHECR 
protons and the CMBR photons (A^ = A* in units of m^/10) . 



Finally, if the out-going pion and nucleon have the same angular momentum, the 
UHECR threshold can be expressed as the following form 

We provide a plot for the dependence of the threshold -E^^™^ on the cosmological 
constant and angular momentums (the in-coming soft photon and out-going nucleon 
and pion) in Fig. 4. 

From the above discussion, one can conclude that a tiny but nonzero cosmo- 
logical constant may provide indeed sufficient corrections to the primary predicted 
threshold j3]. For the observed cosmological constant (which is around the level of 
lO-s^GeV^), if the CMBR possesses a quantum number of the order of 10'^°, the 
threshold will be above the energies of all those observed UHECR particles. The pre- 
dicted threshold should be upgraded to a more reasonable level. Now we can say that 
a possible origin of the cosmic ray threshold anomaly has been achieved. It is the 
cosmological constant that increases the GZK cut-off to a level above the observed 
UHECR events. 
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5 Conclusions and remarks 



In this paper, we have showed a cosmological scenario of constant curvature space- 
time with a homogeneous density perturbation. This is in agreement with the astro- 
nomical observations on supernovae and CMBR that about two thirds of the whole 
energy in the Universe is from dark energy and the lowest order description of the Uni- 
verse may be a de Sitter spacetime. The matter world was dealt with as a perturbation 
around the de Sitter background within the framework of standard cosmology. 

Kinematics in an asymptotically dS spacetime and, in particular, the perturbation 
around the dS background was presented. The exact conservation law of momentum 
and angular momentum in dS is violated by virtue of the additional homogeneous 
density perturbation. Numerical simulations show that, within the range that we are 
interested, the violation of the conservation law is small. A general form of dispersion 
relation for free particles moving in the Universe was obtained. 

The perturbation around dS spacetime has been used to discuss the UHECR thresh- 
old anomaly. We obtained explicitly the corrections of the GZK threshold for the 
UHECR particles interacting with soft photons, which are dependent on the cosmo- 
logical constant. We showed how the threshold varies with a positive cosmological 
constant and additional degrees of freedom of the angular momentums of interacting 
particles. It should be noticed that, for a positive cosmological constant, the theoretic 
threshold tends to be above the energies of all the observed events. Thus, we may 
conclude that the tiny but nonzero cosmological constant is a possible origin of the 
threshold anomaly of the UHECR. 
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